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LOLAS: an optical turbulence profiler in the atmospheric 
boundary layer with extreme altitude-resolution. 
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ABSTRACT 

We report the development and first results of an instrument called Low Layer Scidar 
(LOLAS) which is aimed at the measurement of optical-turbulence profiles in the 
atmospheric boundary layer with high altitude-resolution. The method is based on 
the Generalized Scidar (GS) concept, but unlike the GS instruments which need a 1- 
m or larger telescope, LOLAS is implemented on a dedicated 40-cm telescope, making 
it an independent instrument. The system is designed for widely separated double-star 
targets, which enables the high altitude-resolution. Using a 200 // -separation double- 
star, we have obtained turbulence profiles with unprecedented 12-m resolution. The 
system incorporates necessary novel algorithms for autoguiding, autofocus and image 
stabilisation. The results presented here were obtained at Mauna Kea Observatory. 
They show LOLAS capabilities but cannot be considered as representative of the site. 
A forthcoming paper will be devoted to the site characterisation. The instrument was 
built as part of the Ground Layer Turbulence Monitoring Campaign on Mauna Kea 
for Gemini Observatory. 
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1 INTRODUCTION 

In recent years, a number of efforts have been conducted 
in the domain of astronomical adaptive optics to increase 
the field of view in which the wave front is corrected 
from atmospheric-turbulence-induced perturbations. Partic- 
ularly, Ground Layer Adaptive Optics (GLAO) (|Rigautl 
2002: iTokovininl 12004) is aimed at compensating the wave 
front deformations arising from turbulence very close to the 
ground and not compensating those originating at higher 
altitudes. This idea is based on the facts that the compen- 
sation of lower altitude turbulent layers provides wider cor- 
rected fields of view (|Chunlll998l ) and that turb ulence close 
to the ground is generally the most intense (e.g. jAvila et al.l 
(l2004l )). 

To develop a GLAO system for a given site, it is required 
to have as much knowledge as possible about the optical 
turbulence in the ground layer (i.e. in the first kilometer of 
altitud e above the site). For example, iLe Louarn &; Hubinl 
(2006) recognize that the biggest unknown in their theoret- 
ical analysis of the performance of an infrared GLAO sys- 
tem is the vertical structure of the ground-layer turbulence 
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(C^(h)) and that proper measurements are urgent to vali- 
date their results. 



In view of such necessity, we have developed a method 
called Low Layer SCIDAR (LOLAS), t he concept of which 
was presented by lAvila Chunl $2Q04 ). It is based on the 
Generalised SCIDAR method (GS) but uses a smaller tele- 
scope and more widely separate d double stars as targets . 
The method is explained in § lEgner fc Masciadril (|2Q07h 
reported the use of the GS technique to obtain profiles 
measurements with an altitude resolution of a few tens of 
meters, which was the best GS resolution achieved so far. 
The fundamental difference of that work with the present 
one is that they use the 1.8-m Vatican Advanced Technol- 
ogy Telescope whereas we use a 40-cm dedicated telescope, 
which makes of LOLAS an indep endent monitor. See ^[4] for 
further comments on the work bv lEgner &; Masciadril (|2007l ). 



The experimental setup and data reduction for LOLAS 
is presented in § 03 The first results are shown in § |4] and 
finally, in §[5] we give the conclusion of this work. 
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2 LOLAS CONCEPT 

LOLAS is b ased on the GS concept ([Avila et al.1 Il997l : 
IFuchs et al.l Il998l ). which consists of computing the 
normalised- mean spatial autocorrelation function of short 
exposure-time images of the scintillation pattern produced 
by a double star. The normalisation is performed with re- 
spect to the autocorrelation of the mean image. The result- 
ing normalised autocorrelation map is dimensionally equiv- 
alent to a scintillation index, which is dimensionless. Here- 
after, we will refer to this normalised-mean autocorrelation 
simply as the autocorrelation. 

In the clas s ical SCIDA R (jRocca et al.1 Il974l : 
ICaccia et all Il987l : IVerninl Il992h . the analysis is made 
in the pupil plane of the telescope, which makes it in- 
sensitive to turbulence close to the ground because th e 
scintillation variance is proportional to h 5 ^ 6 (|Roddierlll98lh , 
where h is the altitude of the turbulent layer (acting as 
a phase screen). In the GS the plane of the detector is 
made the conjugate of a virtual plane (analysis plane) 
located at an altitude h gs of the order of a few kilometers 
below the ground, hence h gs < 0. In this case the distance 
relevant for scintillation produced by a turbulent layer at 
an altitude h is | h — h gs | , which makes the tu r bulen ce at 
ground level detectable. Following I Avila et"aD |l997) and 
references therein, the autocorrelation function obtained 
can be written as 
r+oo 

C gB (r) = / dh K (r, h - M C% (h) + N(r), (1) 

J h gs 

where the kernel K (r, h — h gs ) is the theoretical autocorre- 
lation function produced by a single layer at an altitude h 
with a unit C 2 N and N(r) is the measurement noise. When 
h < 0, Cn(/i)=0. For a double star of angular separation 
p, the Kernel consists of three autocorrelation peaks: one 
centered on the autocorrelation origin and the two oth- 
ers separated from the first one by d r = p\h — h gs \ and 
di = — p\h — h gs \. The altitude of the layer is determined 
from the measurement of d r or di and the solution of the 
corresponding equation. 

The altitude resolution is equal to Ad/p, where Ad 
is the minimal measurable difference of the position of 
two autocorrelation peaks. The natural value of Ad is the 
full width at half maximu m L of the auco rrelation peaks: 
L(h) = 0.78 v / A(/i - h gs ) (jPrieur et alJl200lh . where A is the 
wavelength. However, Ad can be shorter than L if the inver- 
sion of Eq. [1] is performed using a method that can achieve 
super-resolution like Maximum Entropy or C LEAN. Both 
meth o ds have been used in GS measurements (|Prieur et al.l 
l200lh . iFriedl (|l995h analised the CLEAN algorithm and its 
implications for super-resolution. Applying his results for 
GS leads to an altitude resolution of 

Ah = iL = ^ 2 VWEM. (2) 

3 p P 

The maximum altitude, /i ma x for which the value 
can be retrieved is set by the altitude at which the projec- 
tions of the pupil along the direction of each star cease to be 
overlapped, as no correlated speckles would lie on the scin- 
tillation images coming from each star. Figure [T] illustrates 
the basic geometrical consideration involved in the determi- 
nation of /imax • Note that /i ma x does not depend on h gs . 
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Figure 1. Schematic for the determination of the maximum alti- 
tude /imax for which the value can be retrieved. The altitude 
of the analysis plane h gs is represented only to make clear that 
this value is not involved in the calculation of h max . 



The maximum altitude is thus given by 

/imax — 5 (3) 
P 

where D is the pupil diameter. 

In addition to the autocorrelation of the scintillation 
images, the GS calculates the mean cross-correlation of im- 
ages taken at time s separated by a known constant dela y 
At. As expl a ined b y lAvila et al.1 (|200lh ; IPrieur et all <|2004h ; 
I Avila etall (|2006h , the cross-correlation leads to the deter- 
mination of the velocity of the turbulent layers and the 
of the turbulence generated inside the telescope dome. 

LOLAS concept consists of the implementation of the 
GS technique on a small dedicated telescope, using a very 
widely separated double star. For example, for h gs = — 2 km, 
h = 0, A = 0.5 /mi, D = 40 cm and star separations of 
180" and 70", Ah equals 19 and 48 m, while h max equals 
458 and 1179 m, respectively. GS uses a larger telescope (at 
least 1-m diameter) and closer double stars, so that the en- 
tire altitude-range with non-negligible values is covered 
(/imax > 30 km). 

The altitude of the analysis plane, h gs was set to -2 km, 
as a result of a compromise between the increase of scintilla- 
tion variance, which is proportional to \h — h gs \ 5 ^ 6 , and the 
reduction of pupil diffraction effects. Indeed, pupil diffrac- 
tion caused by the virtual distance between the pupil and 
the analysis planes provokes that Eq. [T] is only an approx- 
imation. The larger h gs or the smaller the pupil diameter, 
the greater the error in applying Eq. [T] We have performed 
numerical simulations to estimate such effect. A succinct de- 
scription of the simulations is presented in Appendix [Al 

The pixel size projected on the pupil, d p , is set by the 
condition that the smallest scintillation speckles be sampled 
at the Nyquist spatial frequency or better. The typical size 
of those speckles is equal to 1/(0). Taking the same values 
as above for h gs and A, yields L(0) = 2.45 cm. We chose 
dp = 1 cm, which indeed satisfies the Nyquist criterion d p ^ 
L(0)/2. The altitude sampling of the turbulence profiles is 
Sh = dp/ p. Note, from the two last expressions and Eq. [2] 
that the altitude resolution Ah and the altitude sampling 
Sh are related by Sh ^ (3/4) Ah for h = 0. 
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Figure 2. Optical layout. The dimensions unit is millimeters. 
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Figure 3. Schematic view of the instrument. The drawing is not 
on scale. 



3 INSTRUMENT 
3.1 Hardware 

A 40-cm Meade telescope installed on an equatorial mount 
sends light to a pair of 50-mm focal-length achromatic lenses 
that form on the detector an image of a virtual plane located 
approximately 2 km below the pupil. Figure [2] represents a 
scaled optical layout showing the optical beams of two stars 
separated by 200" as they travel from the telescope focal 
plane to the detector. The beam sizes correspond to the 
telescope-pupil diameter. The diameter of each beam at the 
detector plane is equal to 1.30 mm. In some GS systems a 
field lens is used to shift the image of the analysis plane 
either further away from the telescope focus or closer to 
it, whether the lens is diverging or converging, respectively. 
The optical setup for LOLAS was designed so to avoid the 
use of a field lens. The detector is an Electron Multiplying 
Charged Couple Device (EMCCD) manufactured by Andor 
Technology with 512x512 square pixels of 16 /mi side each. 
The makei[j asserts that the effective readout noise is smaller 
than 1 electron r.m.s. and the quantum efficiency is higher 
than 90%. Frames are binned 2x2. The size of each binned 
pixel on the conjugated plane is 9.8 mm. One can chose 
any size for the active sub- array of the detector to be red. 
Figure [3] shows a schematic view of the instrument. The 
telescope holds a motrised rotator used to align the pixel 
lines along the separation of the double star. The rotator 
holds the mechanical support that carries the rest of the in- 
strument. A sliding stage, which enables precise positioning 
along the optical axis, is screwed on the support and car- 
ries the camera. The same equipment, except for the optics, 
is used for Slope Detect ion and Ranging (SLODAR) obser- 
vatio ns. The SLODAR ([Wilson et al.ll2004l : iButterlev et aD 
2006) exploits Shack-Hartmann wavefront sensor measure- 
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Figure 4. Example of a scintillation image acquired with LO- 
LAS. The exposure time was set to 3 ms. The analysis plane is 
located at h gs = —2 km. The object is the double star "STF 
4" of coordinates 0:2000 = Olh 56(2 and 5 2 ooo = +37° 15'. The 
V-magnitudes of the stars are 5.68, 6.10 and their separation is 
199^7. 



ments of the slope of the phase aberration produced by 
the atmospheric turbulence. This technique is implemented 
by replacing LOLAS focusing optics by a collimator and a 
lenslet array. To easily switch between LOLAS and SLO- 
DAR configurations, a manual mechanism is attached to 
the camera and holds the optics of the two experiments. 
The mechanism is designed to ensure the correct position 
of each optical assembly. Figure [3] shows the LOLAS optics 
barrel in place. Into this barrel the LOLAS focusing lenses 
are mounted. The telescope, camera, stage and rotator are 
controlled by a Personal Computer (PC) with two 3-GHz 
Xeon processors and 1-GB memory capacity, running under 
Linux operating system. 



3.2 Data acquisition and processing 

The out-of-focus pupil images produced by each star are 
centered on the detector. To optimize the acquisition and 
processing speed, only a sub-frame of 256x80 binned pixels 
is acquired. An example is shown in Fig. The exposure 
time of each frame ranges from 3 to 10 ms, depending on 
the wind conditions, and is selected by the user. The typical 
number of frames used to obtain one set of auto- and cross- 
correlations is 30000. The frame rate at 3 ms exposure time 
is 112 frames per second. 

In windy conditions, telescope shaking can cause image 
motion which affects the calculation of the mean image. To 
avoid this, every image is recentered prior to co-addition of 
the frames. In addition, telescope tracking errors are cor- 
rected by updating the telescope position every 200 images, 
using the average image position of the most recently ac- 
quired 200-frames packet. This autoguiding keeps the extra- 
pupil images on the active sub-array of the detector. It has 
been seen that the pupil images slowly change in size. This is 
due to a slow shift of the telescope focus, presumably caused 
by the redistribution of the primary mirror load while track- 
ing the stars. An autofocus system has been developed to 
overcome this problem. On every average image calculated 
with 200 frames, the size of the out-of-focus pupils is moni- 
tored and the camera and optics are moved by acting on the 
supporting sliding stage, until the nominal size of 40 binned 
pixels is recovered. Similarly, the alignment of the stars along 
the pixel lines is measured and eventually corrected auto- 
matically by acting on the motorised rotator. Finally, the 
average flux over the pupil is also monitored in each image 
and its variance relative to the mean flux is calculated and 
stored in the header of the FITS file in which the correla- 
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tions are saved. This information is used to qualify the data. 
Indeed, mean flux variations, due to cloud passages or fog 
condensation, are adverse to the retrieval of profiles from 
the autocorrelations. 

The auto- and cross- correlations are calculated as im- 
ages are being acquired. When a set of correlations are saved 
on disk, an IDL program is automatically executed to invert 
Eq. [T] and retrieve C^(h), using a modified CLEAN algo- 
rithm similar to that developed by IPrieur et al.l (|200ll ) for 
GS measurements. 



4 FIRST RESULTS 

The first results of LOLAS as described in §[3] were obtained 
in September 2007 at Mauna Kea Observatory, as part of a 
collaboration between the Universidad Nacional Autonoma 
de Mexico (UN AM), the University of Durham (UD) and the 
University of Hawaii (UH), under a contract with Gemini 
Observatory. The instrument was installed on the Coude 
roof of the UH 2.2-m telescope. 

Figure [5] shows an example of the autocorrelation map 
of scintillation images obtained with the double star "STF 
28" of angular separation of p — 108" 4, which is considered 
a moderate separation for LOLAS. The number of frames 
used for this example was 30000 and the exposure time of 
each frame was 3 ms. The central peak is the result of the 
summation of the central peaks produced by all the tur- 
bulent layers in the atmosphere, even those that are higher 
than the maximum altitude /i ma x reached by the instrument, 
which in this case is equal to 764 m (Eq. [3}. The correlated 
speckles produced by those layers are separated by a longer 
distance than the pupil diameter, thereby not giving rise to 
lateral peaks inside the autocorrelation-map bounds. Only 
the layers below /i ma x form visible lateral peaks from which 
C^(h) is retrieved. The right-hand-side wing of those peaks 
in Fig. [5] are framed by rectangle A. The most intense peak 
inside that frame corresponds to the ground level and is a 
consequence of the turbulence inside and outside the tele- 
scope tube. Much fainter peaks further apart from the map 
center - thus from higher layers - can be seen inside frame A. 
It is worth noting that the color scale for that figure was cho- 
sen such that the fainter peaks are visible, not caring if the 
stronger peaks appear saturated. The autocorrelation signal 
is contained inside a strip which is vertically centered and 
aligned along the horizontal axis, because the detector pixels 
were aligned along the direction of the double-star separa- 
tion. Hence, the values inside frame B correspond only to 
noise. The standard deviation a ac of those values is used 
as a measure of the noise le vel in th e data . In this case, 
cr ac = 4.5 x 10" 4 . Following [Roddierl (|l98lh . the relation 
between C^(0) and a scintillation variance ai is: 



c£(o) 



<ji 



-5/6 



19.12 S h 



(4) 



For the data shown in Fig [5j <5^ = 19 m, hence cr ac corre- 
sponds to a Cn uncertainty of AC^ = 1 X 10~ 16 m _2/3 . 
Correlation values larger than 3cr ac are considered as being 
actual signal. This sets the minimum detectable C§ value 
to CN min = 3 x 10~ 16 m" 2/3 . Note that Ah is defined as the 
resolution along the line- of- sight. For this case, Ah — 31 m, 



0.01 



0.02 



0.03 



0.04 



0.05 



0.06 



0.07 



Figure 5. Example of an autocorrelation obtained with a mod- 
erately separated star (p = 108"4). The scale is indicated by the 
false-color code at the bottom. As explained in §[2] the autocor- 
relations are dimensionless. 
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Figure 6. Examples of profiles. LOLAS profile was obtained 
from the autocorrelation shown in Fig. in 2007 July 3 at 9:51 
UT. The SLODAR profile was obtained the same date at 10:19 
UT. C2 uncertainties for LOLAS and SLODAR are ±1 x 10~ 16 
and ±3 x 10- 16 m" 2 / 3 . 



for h=0. The altitude resolution at ground level in the ver- 
tical direction is Ahcos(z), where z is the zenith angle. In 
this case z — 35.6° which results in a vertical resolution 
of 25.4 m. The profile corresponding to that autocor- 
relation is shown in Fig. [6] The contribution of turbulence 
inside the telescope tube has been removed from the 
value at ground level. To do so, we used the temporal cross- 
correlation of images taken with a temporal interval of 13 ms 
(26 or 39 ms can also be chosen, depending on local wind 
speed), calculated with the same data as that used for the 
autocorrelation of Fig. [5] This cross-correlation is shown in 
Fig. [7] The correlation triplet that is centered on the cross- 
correlation origin is interpreted as a rising from s tatic tur- 
bulence inside the telescope tube (see lAvila et al.l (|200lh for 
details of that method). The corresponding value of the 
is 8.9xl0 -15 m _2/3 , which is equivalent to 0"37 seeing. 

In Fig. [6] is also shown a profile measured with the 
SLODAR, 28 minutes later than the measurement done with 
LOLAS, using the same stellar source. The uncertainty 
is obtained using the bootstrap method. The altitude reso- 
lution achievable with LOLAS is approximately 5 times bet- 
ter than that for SLODAR. This is clearly illustrated by the 
profiles shown in Fig. [6] Around 200 m, LOLAS separates 
two layers, centered at 150 and 210 m, while the SLODAR 
delivers a single peak centered at 190 m, approximately. 
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Figure 7. Cross correlation obtained with the same data as the 
autocorrelation shown in Fig. [5] 
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Similarly, in the 70 m, LOLAS identifies two layers while 
the SLODAR only one. Note that the sum of the val- 
ues delivered by LOLAS in the latter altitude range agrees 
with the SLODAR value at ground level within 10%. Both 
measurements are free of the effect of the turbulence inside 
the telescope tube. The lower altitude resolution of SLO- 
DAR is related to a wider spatial sampling on the analysis 
plane. Here we used an array of 8 x 8 sub-apertures, each 
of diameter 5 cm, which is five times larger than d p . This 
implies that the photon flux per spatial resolution element 
and hence the instantaneous signal-to-noise ratio of LOLAS 
measurements is five times lower than for the SLODARs, so 
that a five times longer period of integration is required to 
measure the turbulence profile to a given accuracy. Hence in 
practice a choice can be made between greater altitude res- 
olution (LOLAS) or greater temporal resolution (SLODAR) 
of the profile. 

To illustrate the highest altitude resolution that has 
so far been reached with LOLAS, Fig. [8] shows a pro- 
file obtained using as target a 199"7-separation double-star. 
The altitude resolution in vertical direction is 11.7 m and 
ACn = 1.6 x 10~ 16 m~ 2/3 . Note the ability for discerning 
a layer centered at 16 m from that at ground layer. Sim- 
ilarly to the profile shown in Fig. [6] the turbulence inside 
the telescope tube has been removed. To make certain that 
the Cn peak at 16 m is not an artifact produced by the 
inversion algorithm, Fig. [9] shows the corresponding auto- 
and cross- correlations. Only an enlargement of the central 
and right-hand-side peaks is shown. In the autocorrelation, 
the right-hand side peak appears slightly elongated in the 
horizontal direction, suggesting that this peak is in fact the 
result of two layers very close to each other. This is con- 
firmed on the cross-correlation map, where one can clearly 
see the contribution of two distinct layers: peaks 1 corre- 
spond to the ground layer and peaks 2 come from a slightly 
higher layer as the separation of the latter peaks is larger 
than that of peaks 1. The altitude resolution might be im- 
proved if the information contai ned in the cross-correla tion 
map was systematically used ([Egner &; Masciadrill2007f ). 



5 CONCLUSIONS 

An instrument has been created for the monitoring of optical 
turbulence profiles in the atmospheric boundary layer with 
an altitude resolution that can reach 11.7 m. It is a stand- 
alone system that needs only an external power supply, mak- 
ing it a suitable device for long-term monitoring at devel- 
oped and undeveloped astronomical sites. The system incor- 
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Figure 8. Example of a turbulence profile with the highest 
altitude-resolution so far obtained with LOLAS. The data was 
taken in 2007 November 17 at 12:09 UT. The central frame shows 
an amplification of the profile in the low-altitude zone. The 
uncertainty is ±1.6 x 10 -16 m -2 / 3 . 
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Figure 9. Similar to Figs. [5] (for a) and [7] (for b), but obtained 
with a widely separated double star (p = 199" 7). 

porates algorithms that automatically preserve the correct 
image position on the detector, focus and camera rotational 
position. Temporal variations of the pupil- averaged flux are 
also calculated automatically and stored in the output files 
as auxiliary data for data reduction. The computation of 
the auto- and cross- correlation is performed in real time 
and the profile is automatically retrieved immediately 
after the correlations are stored on disk. Numerical simula- 
tions of the effect of pupil diffraction show that it reduces 
the values of the autocorrelation by 10%. Nevertheless, this 
alteration is corrected before applying the inversion algo- 
rithm. The turbulence originated inside the telescope tube 
can be removed from the profile on a post-processing 
step. The first results have shown the very high capabili- 
ties of LOLAS in terms of altitude resolution. The shortest 
vertical resolution is only twice that of the instrumented 
ballo ons, which is of the order of 6 m ([Azouit &; Verninl 
2005). One disadvantage of LOLAS compared to the bal- 
loons, SL ODAR, GS, Mult i- A perture Scintillation Sensor 
MASS) (iKornilov et all l iooiih or Single Star Scidar (SSS) 
Habbib et aU 2006) is that until now it does not provide 
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values of the free atmosphere. However, scintillation gener- 
ated at higher layers is indeed detected with LOLAS. This 
information could be used in two possible ways: perform- 
ing a post-detection spatial filtering using MASS-like masks 
or using SSS-like algorithms to retrieve (Rvalues from the 
central peaks of cross-correlations obtained with a set of 
different temporal intervals. SLODAR is the remote-sensing 
technique that has a vertical resolution closer to that of LO- 
LAS. The advantage of SLODAR over LOLAS is that the 
temporal resolution of the C^(h) measurements is five times 
higher. A first comparison of the profiles delivered by LO- 
LAS and SLODAR have shown a fairly good agreement be- 
tween both techniques, although a more thorough compara- 
tive study needs to be done to draw a definitive conclusion. 
The turbulence profiles shown here are not to be taken as 
representative of the site. A proper monitoring campaign at 
Mauna Kea is currently being carried out and a statistical 
analysis of the results will be the subject of a forthcoming 
paper. 
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APPENDIX A: PUPIL DIFFRACTION EFFECT 

The effect of pupil diffraction has been investigated by nu- 
merical simulations of the normalised autocorrelation of 
scintillation images of a double star, on a 40-cm telescope. 
There is no analytical investigation of thi s effect in the lit - 
erature. We used the SCIDAR Simulator (jCruz et al.ll20Q3h 
which is based on the optical turbulence simulator named 
Turbulent 

Plane waves coming from two stars separated by 1 arc- 
minute pass through three turbulent layers located at alti- 
tudes h of 100, 300 and 700 m above the ground. The three 
layers have equal corresponding to a Fried parameter of 
20 cm. The phase of the waves are distorted by the turbulent 
layers following a Kolmogorov spectrum. Between the layers 
and down to the simulated detection plane, wave propaga- 
tion is taken into account by calculating Fresnel diffraction 
of the complex amplitudes. At ground level, a mask is ap- 
plied to the complex amplitudes, which simulates an unaber- 
rated pupil of a 40-cm Meade telescope, including the central 
obscuration. The intensity of the resulting waves are calcu- 
lated on the detection plane, located 2 km below the ground 
(i.e. h gs = —2 km). The autocorrelation of the intensity map 
is calculated and co- added for 1000 statistically independent 
samples. This co-added autocorrelation is normalized by the 

2 http://www.mpia.de/AO/ATMOSPHERE/TurbuLenZ/tlz.html 
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autocorrelation of the mean intensity. The amplitude and 
shape of the correlation peaks obtained for the correspond- 
ing layers are compared to the theoretical ones that would be 
obtained with an infinite pupil-size. Those theoretical values 
are calculated by considering that the correlation amplitude 
is proportional to C^h 5 ^ 6 and the width of the correlation 
peak is proportional to \f\L. The amplitude of the simulated 
correlation peaks are 3.4%, 8.4% and 9.3% smaller than the 
amplitudes of the theoretical peaks, for the layers at 100, 
300 and 700 m, respectively. The measured autocorrelation 
peaks are corrected by a factor obtained by the interpolation 
of those factors for the altitude of the corresponding layer, 
before applying the inversion procedure which is designed 
for an infinite pupil-size. 



